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1. General Remarks 



In the talk given by one of the authors (A. D.) at the AlHENP-95 Conference in Pisa, 
some recent (i.e. since AlHENP-93) developments in two-loop calculations were reviewed. 
The main results discussed in the talk were published |jl|, 0, || and were obtained in 
collaboration with N.I. Ussyukina [|I|], F.A. Berends and V.A. Smirnov 0], and the other 
author of this paper, J.B. T. 0, ||. 

In ref. [|I| (see also in ref. some new results for two-loop three-point diagrams 
with massless internal particles were obtained, including the diagrams with irreducible 
numerators. The paper continues the project started in refs. ^] and provides 
explicit results for the small momentum expansion of two-loop self-energy diagrams for 
all cases when the smallest physical threshold is zero ("zero-threshold expansion"). To 
do this, general results on asymptotic expansions of Feynman diagrams were applied (see 
e.g. in the reviews 0] and references therein). In the paper [H, we have established 
a useful connection between two-loop massive vacuum diagrams and off-shell massless 
triangle diagrams. This connection is valid for all values of the space-time dimension n 
and masses (which are related to the external momenta squared of the triangle diagram 
via ml = pji i = 1, 2, 3). As a corollary, the analytic result for the e-part (e: = (4 — n)/2) 
of two-loop vacuum diagram was obtained in terms of trilogarithms. For the case of equal 
masses, a new transcendental constant (related to log-sine integral Ls3(|7r)) is shown to 
appear. 

Since it is impossible to give a detailed description of all these results in a short 
contribution, and also because the main results have been already published, we decided 
to concentrate ourselves on an important particular problem related to two-loop vacuum 
diagrams, namely the problem of tensor decomposition. 

Tensor reduction of Feynman integrals containing loop momenta with uncontracted 
Lorentz indices in the numerator is very important for various realistic calculations in the 
Standard Model (and beyond). For one-loop diagrams with different numbers of external 
lines, several approaches and algorithms were developed . For two- loop vacuum and self- 
energy diagrams, the problem was considered e.g. in refs. 0, whilst the three-point 
two-loop case is more complicated and requires results for the integrals with irreducible 
numerators 

The problem of finding general algorithms for the tensor decomposition of two-loop 
vacuum diagrams is interesting because it is connected with the calculation of coefficients 



of the small momentum expansion also for the three-point case [jTT[. In the three- 
point case, the problem is more tricky (even for scalar integrals), because we have two 
independent external momenta to contract with. Some relevant formulae for cases when 
the numerator is contracted with one or two external vectors can be found in refs. f^, |12|. 
We also note that the general result for a special case (when two masses are equal and 
the third is zero) was presented in ref. ||10|| . 
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2. Tensor Decomposition 

In this paper, we shall use the following notation: 

/[something] = J J d"p d"g {something} F (j)^ , , (pq) 
and we are interested in expressing the integrals 



(2) 



in terms of scalar integrals. In eq. (|l|), F {p^ , q'^ , (pq)) is an arbitrary scalar function 
depending on Lorentz invariants of the loop momenta p and q. In the special case when 
this function does not depend on (pq) we shall write it as F (p^, g^). Usually, this function 
is a product of propagators, 

(p^-m^) ' (g^-ma) ' ((p- 9)^-7/13) ' or (p^-m^) ' {q'^-ml^ ' , (3) 

but all formulae we are going to discuss are valid for arbitrary scalar functions. 

Some explicit results for the cases when the integral (0) was contracted with the 
external vector (momentum) k were given in ref. (eqs. (B.IO) and (B.9)). They are also 
vahd when the integrands (|]) are replaced by (pq)) (in eq. (B.IO)) or F {jp',q^) 

(in eq. (B.9)), respectively. Moreover, we note that eq. (B.9) can be generalized to the 
case with two external momenta ki and /c2, 



/ ldy<rqF{p',q') [2(M]^M2(M]^M2(pg)] 



N1+N3- even 
N2+N3- even 

A^i! N2\ Nsl 



(^/2)(Ar-^+Ar2)/2 ("-/2)(^^_,.7V3)/2 {''T'/'^) (N2+N3)/2 

{kiy^{kiy^[2{k,k,)y^ (^/2)(iv,+iv.+iV3-,3)/2 



X E 

2ji+j3=Ni 

'2j2+j3=N2 



Jl! J2!j3! ((iV3-j3)/2)! 



X 



/ / d> d-g F (p^ q') (/) W+A^3)/2 (^2)(A'.+A^3)/2 



where {a)j = r(a + j)/r(a) denotes the Pochhammer symbol. If {Ni + A^3) or (A^2 + ^3) 
is odd, the integral on the l.h.s. is equal to zero. 

Since we have no external momenta in (|^), only tensor structures constructed of metric 
tensors may be involved in the tensor decomposition. Therefore, (0) vanishes if = 
A^i + N2 is odd. For even N, the suitable tensor structures should be symmetric with 
respect to two subsets of indices, (/xi, . . . , /zatJ and (ai, . . . , a^^). All independent tensor 
structures can be constructed in the following way. Let us take a product of ji metric 
tensors (7^,^^. , j2 tensors Qa^ak and js tensors fi'^io-j. (we should remember that the conditions 
2ji + js = Ni and 2j2 + js = N2 should be satisfied). Now, let us consider all possible 
permutations of /i's and all possible permutations of cr's producing distinct products of 
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metric tensors, and take the sum of all these terms (with the coefficients equal to one). 
So, the obtained tensor structure is symmetrized in /x's and cr's, and it is 



+permutations. 

For given A^i and N2, the number of independent tensor structures (1) is 



+ 1, 



where 



IN- 
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T(iVi,iV2)=min([iiVi],[iiV2^ 

is the integer part of and the number of terms on the r.h.s. of 

Ni\ N2\ 



3l]2]3 



(5) 



(6) 



IS 



(7) 



Note that due to the conditions 2ji + J3 = Ni and 2j2 + js = A^2, at given Ni and N2 the 
tensor structures (H) are completely defined by one index, J3. 

So, the result for the integral (0) (for even Ni + N2) should look like 



^N2 



2ji+j3=Ni 

2j2+j3=N2 



where /jijajs ^"^^ some scalar integrals which are to be found. A standard way to define 



^J1J2J3 



is to consider all independent contractions of (||) (e.g. to contract it with each of 
the structures (^). In such a way, one needs to solve a system of T{Ni, N2) (see eq. (|^)) 
linear equations. 

Due to power-counting reasons, the scalar integrals Ijijzja should be linear combina- 
tions of the integrals (|I]) with scalar numerators, carrying the same total powers of the 
loop momenta p and q. Therefore, we can re-write eq. (|^) as 



P^il ■ ■ ■ P/iJV]^ QcTl ■ ■ ■ Qc 



■JV2 



= E {jl,j2,j3} E (^jlj2j3j[j'2j'3 Cjif.ji, I 

2ji+j3=Ni 2j[+j!^=Ni 

2j2+j3=N2 2j^+j'^=N2 

where Cj^j2i3 defined in (|^, while ^ju^iayjiais '^^^ considered as the elements of the 
"decomposition" matrix (T x T, see eq. (|])). Contracting eq. (^ with all possible tensors 
(H), we get a "column" of Cj^j^j.J [{p^Y^ {q^)^^ {pq)^^] on the l.h.s. On the r.h.s., we get the 
"contraction" matrix x with the elements defined as 

Xju2j3;iij^j^ = contraction ({ji, j2, js} , {j'id'2^ J3}) • (10) 

When speaking about matrices, we understand that only one of j-indices is relevant in 
each set (ji, 72,73) (e.g., 73). For given Ni and the values of js can be either even (if 



(9) 
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A^^i and N2 are even) or odd (if Ni and N2 are odd). In order to use the matrix notation 
in ordinary form, we introduce the generahzed index j = (ji, 72,73) (and j' = U'li j'21 f'i) 1 
etc.) which takes values from 1 to T{Ni, N2). In all "non-matrix" formulae, however, 
we would prefer to keep all the indices ji, 72, js explicitly. So, in the matrix notation, a 
corollary of eqs. (pD-(p!OD is 

J2 Xjj" <Pj"j' = Sjj' ■ (11) 
j" 

In other words, the decomposition matrix is nothing but the inverse contraction matrix, 
X~^- Since Xjj' is symmetric (see (pUj)), is symmetric as well. 

Thus, the problem is how to find all 4'jij2j3;j[j!2j'^ given A^i and N2. One way is to 
use recurrence relations. We note that the following simple formulae of contracting tensor 
structures (^ with respect to two indices can be derived: 

<xiV2 { Ji ' J2 J3} = {n + 2ji + 2j2 + J3 - 1 ) { Ji , J2 , Js - 1 } + (js+l) { ji - 1 , J2 - 1 , Js + 1 } ,( 12) 

{ji, J2, Js} = (n + 2ji + 2j3-2) j2, j3} + 2 (J2 + I) {ji, J2 + I, J3-2} (13) 

(and also an analogous formula for contraction with S'o-jvj-io-jvj)- Using these contractions, 
the following recurrence relations for ^jijajayjiais obtained: 

A^iA^2 {{n + Ni + N2- J3 - 1) (Pnj,j,-j[j!^f^ + (js - 1) (pj,j,j,-j[+i,j'^+ij'^-2} 

= h 0iij2,i3-l;i;,i^4-l ' (14) 

Ni (iVi - 1) {{n + Ni+j!^ - 2) ^ju^jauli^i^ + 2j2 (l)j,j,j,-j[-i,j^-ij'^+2} 

= 2 ji 0ji-lj2j3;j;-lj^,j^ , (15) 

and also a relation similar to (0) but with interchanged indices, 1^2. In eqs. (|T^-(0), 
it is understood that, whenever any of j's becomes negative, the corresponding tensors 
J2, js} and 0's should be taken equal to zero. 

If we increase A^i and/or N2 without changing the number of tensor structures (|^), 
it is enough to have only relations (p!^)-([T5|) to express "higher" 0's in terms of "lower" 



0's. If on this step one extra tensor structure appears (for example, when we go from 
odd A^^i and N2 to {Ni + 1) and (A^2 + 1)), we need one more relation between </)'s. This 
extra relation can be obtained by contracting (^) with k^^^ . . . k^^^ k""^ . . . fc'^^a and using 
eq. (B.IO) of ref. g, 

S <Pnnh^3[3',,i, = 2W+7V,)/2 {^n/2)<N,+N,)/2 ' ^^^^ 

The values of the lowest 0's can be taken from the results at A^2 = and A^2 = 1, 



iVi 



^ ^^,0,0 ^ / {p^r^'^ , (17) 



2^1/2 (n/2)^,/2 I 2 
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Ni - 1 



0,1 / 



ip^Y^^'^y^pg)] . (18) 



2W+i)/2 (n/2)(;v,+i)/2 I 2 
Using the recurrence procedure, it is also possible to obtain some less trivial formulae, for 
N2 = 2 and A^2 = 3 and arbitrary A^i (such that A^i > iV2): 



2(A.,+2)/2 (n/2)(^,+2)/2 {n - 1) 



X 



,1,0 



(n + N,- I)(p2)^i/2g2 _ J^^^p2YN^-2)/2^p^y 



'Ni-2 
+ ^^7^,0,2 



n(p2)(^i-2)/2(p^)2_^p2)JV,/2^2 



(19) 



Pt^l 



X 



■ PmaTi ^o"2 ^crs 

Ni-l 



2 

+ 



1,1 / 



2(^i+3)/2 (n/2)(^,+3)/2 (r^-1) 

(^ + iVi-2)(p2)(^-l)/2g2(p^) _ (^^ 



l)(p^ 



{pqf 



,0,3} I Un + 2)(p2) W-3)/2(p^)3 _ 3(p2)(7Vi-l)/2^2^p^) 



(20) 



'^Jlj2j3;Jii2J3 



is constructed which provides a 
It is difficult, however, to generalize 



Thus, a recursive procedure of calculating 
general algorithm to calculate the tensor integrals 
eqs. (p!7|)-(|20D to the case of arbitrary N2. 

There is, however, another approach to this problem which does not involve using 
recurrence relations. Let us look again at the decomposition (|^), remembering also how 
the integrals / are defined (p. Then, let us contract eq. (||) with . . . fcf^^fcg^ . . . k'^'^^ , 
and multiply it by (^1^2)^^ ^i)- Now, we can integrate the resulting expression over 

ki and /c2 (not over p and g!), using the formula (^. As a result of this integration at 
different values of N^, we get T{Ni,N2) independent equations. Then, considering the 
scalar integrals I [{p'^y^ig'^y^ipq)-'^] as a basis, we get a set of relations for 0's. Finally, 
taking some linear combinations of these relations, we get a Kronecker symbol on the 
r.h.s. This means that we have found the elements of the matrix inverse to (pjj' which is 
nothing but the contraction matrix Xjj' (see eqs. (p!OD-(pJ])). The obtained expression is 



Xj 



lJ2J3;JiJ2J3 



jV-32ij3' 



E 

/=0 



N,\ m i-iy (J3+J3- 2/)! ((n + N, + N2-2l + 2)/2)i 



l\ (j^-2/)! ((is + - 20/2)! (n/2)(,3+,.^20/2 

x(^/2)(iVi+j--2/)/2 (^/2)(JV2+i^-2/)/2 • (21) 

Thus, we have obtained an explicit general expression for the inverse of the decompo- 



sition matrix, {(p For the cases of interest, the matrix ( PT] ) can be inverted by use of 

computer systems for analytical calculations. Although it would be nicer to get an analo- 
gous result for the "direct" matrix cpjj', the presented approach can also be considered as 
a general solution to the problem of tensor decomposition of two-loop vacuum diagrams. 
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